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Abstract. Future wired and wireless communication systems will employ pure or combined Code Division 
Multiple Access (CDMA) technique, such as in the European 3G mobile UMTS or Power Line Telecommu- 
nication system, but also several 4G proposal includes e.g. multi carrier (MC) CDMA. Former examinations 
carried out the drawbacks of single user detectors (SUD), which are widely employed in narrowband IS-95 
CDMA systems, and forced to develop suitable multiuser detection schemes to increase the efficiency against 
interference. However, at this moment there are only suboptimal solutions available because of the rather 
high complexity of optimal detectors. One of the possible receiver technologies can be the quantum assisted 
' computing devices which allows high level parallelism in computation. The first commercial devices are 

estimated for the next years, which meets the advert of 3G and 4G systems. In this paper we analyze the 
error probability and give tight bounds in a static and dynamically changing environment for a novel quan- 
tum computation based Quantum Multiuser detection (QMUD) algorithm, employing quantum counting 
algorithm, which provides optimal solution. 
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1. Introduction 

> 

The subscribers of next generation wireless systems will communicate simultaneously, sharing the same 
frequency band. All around the world 3G mobile systems apply Direct Sequence - Code Division Multiple 
Access (DS-CDMA) promising high capacity and inherent resistance to interference, hence it comes into 
the limelight in many communication systems. Nevertheless due to the frequency selective property of the 
, channel, in case of CDMA communication the orthogonality between user codes at the receiver is lost, which 
leads to performance degradation. Single-User detectors were overtaxed and shown rather poor performance 
even in multi-path environment [Q. To overcome this problem, recent years multiuser Detection (MUD) 
has received considerable attention and become one of the most important signal processing task in wireless 
communication. 

Verdu JjJ has proven that the optimal solution is consistent with the optimization of a quadratic function, 
which yields in MLSE (Maximum-Likelihood Sequence Estimation) receiver. However, to find the optimum 
is an A/P-hard problem as the number of users grows. Many authors proposed sub-optimal linear and non- 
linear solutions such as Decorrelating Detector, MMSE (Minimum Mean Square Error) detector, Multistage 
Detector, Hoppfield neural network or Stochastic Hoppfield neural network jy, ||, ||, and the references 
therein. One can find a comparison of the performance of the above mentioned algorithms in [Q. 

Nonlinear sub-optimal solutions provide quite good performance, however, only asymptotically. Quantum 
computation based algorithms seem to be able to fill this long-felt gap. Beside the classical description, which 
we recently use, researchers in the early 20 th century raised the idea of quantum theory, which nowadays 
becomes remarkable in coding theory, information theory and for signal processing Q . 

Nowadays, every scientist applies classical computation, using sequential computers. Taking into account 
that Moore's law can not be held for the next ten years because silicon chip transistors reach atomic scale, 
therefore new technology is required. Intel, IBM, AT&T and other companies invest large amount of research 
to develop devices based on quantum principle. Successful experiments share up that within 3-4 years 
quantum computation (QC) assisted devices will be available on the market as enabling technology for 3G 
and 4G systems || § . 

This paper is organized as follows: in Section |^ we shortly review the necessity of multiuser detection, as 
well as the applied quantum computation method is shown. In Section ^ the proposed quantum multiuser 
detector model is introduced. Furthermore, in Section we give and proof a tighter probability error for 
our model. Finally we conclude our paper in Section |^. 

2. Theoretical Backgrounds 

2.1. Multiuser Detection. One of the major attributes of CDMA systems is the multiple usage of common 
frequency band and the same time slot. Despite the interference caused by the multiple access property, 
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the users can be distinguished by their codes. Let us investigate an uplink DS-DCDMA system, where the 
i th symbol of the k th (k = 1, 2, . . . , K) user is denoted by b k (i), b k (i) S {+1, -1}. In DS-CDMA systems 
an information bearing bit is encoded by means of a user specific code with length of the processing gain 
(PG)Q. In case of uplink communication we assume perfect power control. In the receiver side it is not 
required synchronization between input signals and user specific codes, however we make our decision on 
symbols. Applying BPSK modulation, the output signal of the k th user, denoted by qk(t), is given as 

qk{t) = \/Ekb k (i)s k (t), 

where E k and s k (t) are the energy associated to the k th and the user continuous signature waveform, 
respectively 

PG-l 

Sk(t) = ^2 c kti)9k{t - jT c ), 

3=0 

T c denotes the time duration of one chip, c k (j) is the j th chip of the code word of subscriber k and g k (t) refers 
to the chip elementary waveform. We investigate a one path uplink wideband CDMA propagation channel. 
The channel distortion for the k th user is modelled by a simple impulse response function h k (t) = a k 8{b — T k ), 
where a k and r k are the path gain and the delay of the k th user, respectively [§). They are assumed to be 
constant during a symbol period of T s . This model contains almost all elements of a typical WCDMA 
channel except multipath propagation, which was omitted to simplify the explanation of the new quantum 
computation based multiuser detection scheme. However, based on the results of the present paper, multipath 
propagation can be included into the channel model on a very simple way. 

The received signal is the sum of arriving signals plus a Gaussian noise component and thus can be written 
as follows: 

K K 

r(t) = ^2 h k (t) * q k (t) + n(t) = ^ ^fE^a k b k s k {t - r k ) + n(t), 

k=l k=l 

where K is the number of users using the same band, n(t) is a white Gaussian noise with constant N$ 
spectral density. 

Unfortunately, the search for the globally optimal MUD p| usually proves to be rather tiresome, which 
prevents real time detection (its complexity by exhaustive search is 0(2 K )). Therefore, our objective is 
to develop new, powerful detection technique, which paves the way toward real time MUD even in highly 
loaded system. Since classical multiuser detection schemes only try to minimize the probability of error in 
noisy and high interference environment, they, even also optimal solutions, can commit an error. Actually, 
these classical approaches make compromize between computational complexity, probability of error and 
time barrier required for efficient working. On the other hand, QMUD provides for typical CDMA systems 
BERw 1CU 10 and it is able to indicate undistinguishable decision situations for correction by higher layer 
protocols ||. 

2.2. Quantum Computation Theory. Quantum theory is a mathematical model of a physical system. 
To describe such a model one needs to specify the representation of the system. According to the axioms 
of quantum mechanics, every physical system can be characterized by means of its states \<p)q in the Hilbert 
vector space over the complex numbers C, whereas a physical quantity can be described as Hermitian operator 
U = U\ respectively. 

In the classical information theory the smallest information conveying unit is the bit. The counterpart 
unit in quantum information is called the "quantum bit", the qubit. Its state can be described by means 
of the state \<p), \tp) = a\0) + /3|1), where a, (3 € C refers to the complex probability amplitudes, where 
\a\ 2 + |/3| 2 = 1 JlC], ||. The expression \a\ 2 denotes the probability that after measuring the qubit it can 
be found in computational base |0), and |/3| 2 shows the probability to be in computational base |1). In 
more general description an AT-bit "quantum register" (qregister) \(p) is set up from qubits spanned by |i) 
i = . . . (N — 1) computational bases, where N = 2" states can be stored in the qregisters at the same time 
|J, describing 

N-l 

(2.1) \<p) = <Pi\i) <Pi G C, 

where N denotes the number of states and Vi ^ j, = 0, = l,J2 Wi\ 2 = 1> respectively. It is worth 
mentioning, that a transformation U on a qregister is executed parallel on all N stored states, which is called 
quantum parallelization. 



Say ket tp (using Dirac's notation). 
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3. Quantum multiuser Detector Employing Quantum Search 

For the optimal decision it would be necessary a fully comprehensive knowledge about the symbols sent 
by all the subscribers in the coverage area of base station, a realization of the delay and noise, which is 
typical for a particular communication channel and all the user specific codes. All this information cannot 
be stored in a single database, which should be build just one time at all. In pursuance of detection one 
compares the quantized received signal with the content of the database. This task can be done efficiently, 
employing the Grover |yl [l^] database search algorithm (£/), which is proved to be optimal |l3| in 0(sfN) 
steps, where N denotes the length of the database. 

The optimal number of iterations is depending on the initial angle 9 between the initial state I7) and |a), 
as well as the number of the identical entries M in the database. In it was shown that the rotation of 
state I7) to the desired state after d evaluations of Grover operator (Q) is 

(3-D g%)=c ^i±i e y a)+sm (ii±i, 

whereas in fl3| it was proved that the optimal number of iterations I is given as 
(3.2) d = floor 



arccos a / 77- 



d—x.h — >d—x 

In case of M <C N or even it exists no multiple entries in the qregister the angle 
(3.3) 

where the optimal number of evaluations is equal to d ~ f ^/ ^ |l3|| ■ It is worth mentioning that in case of 
N = AM, the initial angle 8 = 60°, which leads to only one evaluation. 

In case of quantum detection, however, it is not required to know the exact position of the representative 
of the received signal in the qregister, only whether is it contained in the database or not. For this purpose, 
the so called Quantum Counting based on Quantum Fourier Transformation (QFT) designed for phase 
estimation is well suited. In the following part of this section we introduce the buildup of the applied 
databases, qregisters. Furthermore, we shortly discuss the quantum counting problem and our way to deal 
with quantum multiuser detection problem. 

3.1. Applied Databases and Qregisters. Let us assume we would like to detect the k th user's symbol, 
denoted by To design a suitable multiuser detector all of the user specific code words have to be 
known at the base station, collected in a (PG x K) sized database C = [c%, C2, . . . , ck]- In addition we create 
a qregister |<^), which contains all the possible received signal configurations \x) without bk 

JV-l 

(3.4) \<p) - J2 <Pi(x)\x), 

i=0 

where \x) are the computational base states. Since every possible noise and delay states for a given trans- 
mitted bit is involved in the prepared register, it is superfluous to send any phase information, which allows 
non-coherent detection. The first {K — 1) qubits represent all the combination of the symbols transmitted 
by all the interferer, e.g. if the value of (K — j) th qubit in a given \x) is equal to 1, it shows that the user 
[K — j) has been transmitted a symbol denoted by bx-j — 1, otherwise if 6fc = —1 the qregister has an 
entry 0. For a proper utilization of the qregister the next two parts in \x) is reserved for all quantized state 
of the noise n(t) and the delay r: < r < T s . For the sake of simplicity we assume to apply a linear 
quantization n = i ■ ur, where i G [— N n ,+N n ] and ur refers to the quantization step. It is remarkable, 
if the probability density function (pdf) of the target quantity is known, such as by the Gaussian noise, 
an adequate nonlinear quantization could be used, whereby most significant values can be sampled denser, 
which can decrease the number of required qbits for description of noise. The values of the delay t can be 



generated in a straitforward manner, r = j ■ tr, where j could be j € 
description the size of qregister |<p) is given by 



With the above mentioned 



(3.5) N qreg = [IdiV] =(K-1 



\d(2N n + 1) + Id I — + 1 



where [•] denotes the round up operation to the next integer. We consider two independent qregisters 

and f° r every user k according to sent information symbol +1 or — 1. From this point we just refer to 

|<£>;), where j = { — 1, 1}, however, the operations on each qregister can be done in the same manner. Thus, 
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\<Pj) are the databases as well as detection means to check whether the received signal can be found in the 
databases [Q. 

3.2. Quantum Counting. Performing multiuser detection it is obvious not to look for a pattern in a 
database itself but to indicate whether or not a solution even exists, consequently, the number of solutions 
have to be summed up. This quantum counting algorithm can be accomplished with combining the Grover 
quantum search algorithm mentioned in Section with the phase estimation based on the quantum Fourier 
transformation [ p^[ . The quantum counting algorithm is able to estimate the eigenvalues of a unitary 
operation U, or in this case the Grover iterations (Q) 

(3.6) 0| 7 )=e^|7), 9 = 2*5, 

which is in a close connection with the demanded number of solutions M on the search. In addition to the 
original Grover circuit a new register \t) of I new qubits should be feed to the quantum counting algorithm, 
whose state is initialized to an equal superposition of all possible inputs by a Hadamard gate, as in Figure [|. 
The vector (7), which contains all the possible signal states with equally distributed probability amplitudes, 
is build up in the same way. The original Grover block Q becomes a controlled Grover block, where the 
transformation (a rotation with angle 9, such as in the original Grover algorithm) will be only evaluated if 
the control signal equals to one. The goal of quantum counting is to estimate 6 to m bits of accuracy, with 
probability of success 1 — e. The state of the register \t) after the controlled Grover iterations is 

oi-l 
k=0 

which is exactly the QFT belonging to 9. Hence, to estimate 6 an inverse quantum Fourier transformation 
is executed on ( |3.7| ) before the measurement. In case of S > after a measurement, the demanded pattern 
is contained in the qregister \<Pj), and if it is equal to zero, it is not included in. 



4. The required number of additional qubits 
The required number of qubits I 

(4.1) 

in \t), and hereby the required number of Grover blocks Q can be calculated in function of desired accuracy 
2~ m of 5 and the probability e of false detection. This error e is occurred if 5 is not in form of z/2, where 
z € 0, 1, . . . ,2 l . Hence it can be regarded as the probability of false detection. 









I = m + 




+ ld7T 









Lemma 4.1. The accuracy m of 5 can be given as 
(4.2) sup{m} 



V 



qreg 



1 



Proof. Fortunately, we are able the upperbound m in (4.1). From the accuracy point of view the worst case 
is occurred if M = 1, because the initial angle 9, yielding form (3.3), is the smallest. Hence, m should be 
chosen as large as possible, which leads to 



= 2 



= 2 



1 



2 N vr 



= 2\ 



sup{m} 



Y 



1 



□ 

Therefore, we have a direct connection between the detection error probability and the required length of 
register \t). However, before the measurement an inverse discrete QFT is done, which adds additional states 
beside \<p) with a given probability, if 9 is not an integer power of two, that could lead to false detection 
in the receiver. Hence, the expression p — [id (2 + ^)] |l4j in ( |4.2| ) is just a approximation. In general, 
p is needed to represent 9 more precisely, and thus to decrease the probability of false detection, but the 
accuracy of the estimated angle 9 is 2~ m in the future too. This means, that p only influences the probability 
amplitudes, which leads to be enough measuring the first m most significant bits! 

4.1. Methods to improve the measurement. To improve the proper error probability of detection, there 
are some possibilities. 

1. One should chose 9 — 2 L , where L £ {N}. Unfortunately, except the trivial case, where 9 = 0, which 
means the desired state \r k ) is not in the qregister, it happens very rarely. 
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FIGURE 1. The quantum counter circuit 



2. Our goal is to distinguish the case 9 = from 9 0. It is known from the previous point, that the 
detection of 9 = is certain. Consequently, if 9 ^ a false detection is occurred if and only if we 
decide for |0) after the measurement. For a more precise estimation of the error probability P errori one 
should extend the size of the qregister instead of p bits, with C more bits. We will show that based on 
our results C can be much smaller than p at a given P err or , where 



(4.3) 



= P(6 = 0\9 ^ 0) • P(9 ^ 0) + P(9 ^ 0\9 = 0) • P(9 = 0), 



where 9 refers to the estimated phase. The second part of the right side in (4.3) is equal zero, because of 



the previous analysis. Assuming an unknown a-priori probability (O) becomes simpler 

Perror < P(6 = O\0 ^ 0). 

The state of the examined \t) qregister before measurement is given as 



^ 2'-l 

Wi E 



k=0 



/ol ^ 



i=0 



IQFT corresponding to \k) 



^EE^ 



(4.4) =^fEE e ^ fc(5 ^ ) N) ^EE^'^w. 



k i i k 

where < 5 < 1 is a real number. Furthermore, the probability amplitude of the state \i) is 

2 ! -l 

a,; 



(1) 1 ^ e ^-*)l E (^-*)) k , 



which is a geometrical series 
(4.5) 
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The demanded probability of error can be described as the sum of over all probability amplitudes given in 

2 C -1 

(4.6) P error < P(9 = 0\9 £ 0) = l««(0| 2 , 



i=0 



form which the parameter C is calculable at a given P e , 



4.2. Upper bound of error probability. While C can not be expressed explicitly from (4.6), hence, in 
case of adaptive receiver structure it would be very useful to give less complex way of determine C . For this 
purpose an upper bound has to be given for (4.5). It is considerable that with increasing the size of the 
qregister with C new bits, the sum in (4.6) may also greater, however, due to more accurate representation 
of 9 the probability amplitude belonging to the false vectors diminishes better and better. 



4.2.1. Numerator of ■ For the numerator of (4.5) in p4| was shown an unambiguous upper bound 
e j27r(<52 -l) 2 ; which is simple the diagonal of unit circle \&* a \. We show now, that under certain 
conditions a tighter bound can be given. 

Theorem 4.2. For a fixed 5 the numerator of the sum of the probability amplitudes is 

max {tt (2 m+c+1 S - noor(2 m+c+1 (5)) ,n-n (2 m+c+1 S - floor(2 m+c+1 )) } . 
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Proof. Unfortunately, the numerator of ( |4.5| ) is not monoton on 5 and I, hence, the maximum has to be 
searched which is less than 2. The unit circle is divided into two regions, where < |a|mod27r < tt and 
tt < |a|mod27r < 27r, where we should find 

max{|a|mod7r} = max{a(7, 5)}, 

1,6 1,6 



max{7r — |a|mod7r} = max{7r — a(l, 5)}, 

1,6 1,6 



respectively, whereby we mapped |a| in the range [0,7r], where \a\ strictly monoton. To get a proper upper 
bound the maximum of the two regions for a fixed 5 



(4.7) 



[(2tt 1 2*5 - l\) mod] tt = [(2 \2*5 - l\) modi] tt. 



The case of, when 2*6 - I < in ([D]), yields in I > 2*5, which can not come true, inasmuch as we sum up 
to 2 C — 1, which is always less than 2*. Hence, the absolute value operator can be neglected. The opposite 
case, where 2*5 — I > in (4.7) the non-integer part becomes tt (2* +1 5) modi because of 21 is an integer. 
These two considerations lead to 



(4.8) 



{tt (2 m+c+1 5 - floor(2 m+c+1 £)) , tt - tt (2 m+c+1 5 - floor(2 m+c+1 )) } . 



4.2.2. Denominator of (4-t). Also in |14j was shown, that the denominator of (4.5) 
always greater than in the range [— 7T, tt]. 



□ 



= le^l is 



Theorem 4.3. Employing the Grover search algorithm in quantum counting the denominator of (4-5) be- 

2V2\/3\ 



Proof. For the proof, we should first evaluate /? m ax, which is 



A. 



n / £■ "mm 

= mmmax27r i) may — 

16 V 2* 



Form (^^), where the maximal value of 9 = ^, resulting from employing the Grover searching in the quantum 
counting algorithm, i5 m ax = \- Furthermore, it is obvious, that ^ m i n = 0, which two considerations lead to 



2n[- 4 -0 



2' 



(i c - 
Omin J* 

where J m ; n is very close to zero and the second part greater than zero except for C > t — 2, which is not 
reali zed almost surely. Following the above mentioned ideas, \(3\ < 5. Applying /3 max the denominator of 
(4.5) becomes 



(4.9) 



2V2 



2,{5- ¥ 



= 4V2(5-^ 



□ 



4.3. Simulation of new bounds. To verify our results we performed a computer simulation of our new 
probability of error upper bounds and compared them with them of the model given in ]l4| , which is denoted 
by Model 1 in Fig. |^. For a static environment we propose Method 3, where an exact calculation of ( |4.6| ) 
is done, whereas Method 2 refers to our upper bound for dynamically changing system, following Lemma 1, 
Theorem 4-2 and Theorem 4-3. In Fig. ||the probability of errors of the outcomes of the three methods are 
depicted against the number of additional bits C. It is noticeable a remarkable difference between Model 
1 and Model 2 & 3 especially at low values of C . Conversely, the likewise less computational complexity 
loaded Model 3 shows a quiet close result to the theoretical one Model 2. It is also observable that the error 
vanish with very few new bits in the qregister. 
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FIGURE 2. The probability of error vs. additional bits 
5. Conclusions 

In this paper we presented a quantum computation based multiuser detection algorithm, which involves 
a modified quantum counting algorithm, employing the Grover's quantum search method. We showed 
an exact computation method for the error probability, which can be obtained in a static environment. 
Furthermore, we also gave an fast computable approximation of the upper bound of error probability very 
close to the theoretical one, which is usable in adaptive receiver. The new method utilizes one of the possible 
future receiver technologies of 3G and 4G mobile systems, the so called quantum assisted computing. QMUD 
provides optimal detection in finite time and complexity when classical methods can achieve only suboptimal 
solutions. The proposed algorithm has strong resistance against MAI and noise. 
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